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The transport coefficients of quark gluon plasma are calculated on a lattice 16 3 x 8, with the pure gauge models. 
Matsubara Green's functions of energy momentum tensors have very large fluctuations and about a few million 
MC sweeps are needed to reduce the errors reasonably small in the case of the standard action. They are much 
suppressed if Iwasaki's improved action is employed. Preliminary results show that the transport coefficients 
roughly depend on the coupling constant as a~ 3 (g). 



1. Introduction and Formulation 

Quark gluon plasma(QGP) is expected to be 
realized in high energy heavy ion collisions in near 
future and thought to play an important roll at 
the very early stage of the universe. Phenomeno- 
logically it will be treated as a fluid and it will be 
very useful to calculate the transport coefficients 
from the fundamental theory of QCD. 

The calculation of the transport coefficients of 
QGP has been done within the Kubo's linear re- 
sponse theory According to this theory, the 
transport coefficients are calculated by the space 
time integral of retarded Green's function at finite 
temperature. It is generally difficult to calculate 
the retarded Green's function at finite tempera- 
ture directly. Then the shortcut is that the Mat- 
subara Green's function is calculated and by ana- 
lytic continuation retarded Green's function is ob- 
taied. For the analytic continuation the spectral 
representation of the both Green's functions is 
used, for which the following ansatz is assumedpl, 



r), 



7T (m — lu) 2 + 7 2 (to + lu) 2 + 7 : 

where 7 is related to the imaginary part of self 
energy. From these parameters, the transport co- 
efficients are calculated as, 



a(A,j,m) = 2A 



27TO 



(7 2 + TO 2 ) 2 



here a represents heat conductivity, bulk and 
shear viscosities. 

We notice that if m = or 7 = 0, transport 
coefficient becomes zero. In order to fix these pa- 
rameters, at least three independent data points 
in temperature direction are required for the Mat- 
subara Green's function. 

In the pioneering work of Karsch and Wyldpf, 
they performed lattice QCD calculations on 8 3 x 
4 lattice and the resolution was not enough to 
determine these three parameters independently. 

2. Matsubara Green's Function 

We carry out the simulation on 16 3 x 8 lattice 
for U(1),SU(2) and SU(3) pure gauge theories. In 
the gauge theory, the energy momentum tensor is 
expressed as follows, 



T M „ = 2Tr{F ll<7 F v 
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where F^ are field strength tensor defined by pla- 
quette variables as U Xjflu = exp (ia 2 g F^ v ) . From 
the plaquette variable F^ is obtained either by 
taking log or by expanding with respect to ga 2 . 
We call the former method as 'diagonal method' 
and a latter one as 'perturbative method'. 

We have started the simulation from compact 
U(l) case, because it has two phase, the confine- 
ment and the Coulomb phase which are separat- 
eded at (3 ~ 1.0, and is a good place to test 
the formulation. In the following we denote the 
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Matsubara Green's function of energy momentum 
tensor T^ v as G^(T), 

In Fig. 1, we show how the error decreases 
with number of MC sweeps in confined and de- 
confined phase for U(l) case. In this work the 
mesurements are done with every sweep. The dif- 
ference in the fluctuations of G M „ (T) between the 
confined and deconfined phase are observed. In 
the confined phase, we could not reduce the ratio 
Errors/ < G(T) > less than unity within more 
than million data even at T — 2. Similar situa- 
tion is observed in SU(2) case. We interpret this 



In order to decrease the fluctuation, we try to 
use an improved action proposed by IwasakiQ for 
the SU(3) case. The phase transition region for 
the Iwasaki's improved action on 16 3 x 8 lattice 
is ~ 2.7 — 2.9 and we started our simulation 
at = 3.3 . In Fig. 2 we have also shown how 
errors decrease with number of MC sweeps for the 
improved action. The fluctuation is much reduced 
by the use of the improved action. From Fig. 2 we 
expect that within a few hundred thousand data, 
we could make the ratio Error(T)/ < G(T) >< 1 
and get the transport coefficients of QGP with 
reasonable accuracy. 



Figure 1. Error as a function of number of MC 
sweeps. 



as follows; the energy momentum tensor in the 
confined phase should be written by the hadron 
fields. Namely used in this calculations are 
not good operators in the confined phase. There- 
fore in the following we calculate the transport 
coefficients only in the deconfined phase. 

The fluctuations of G^{T) for U(1),SU(2) and 
SU(3) cases in the deconfined phase are shown 
in Fig. 2., where transition point is ~ 6.05 for 
SU(3) U and cross over region is ~ 2.4 for SU(2) 
for Nt — 8. The fluctuation becomes larger as we 
proceed to SU(2) and SU(3). And as the fluctua- 
tion increases as T becomes large, we need about 
million data for SU(2) = 3.0 case to make the 
ration Error (T)/ < G(T) >< 1 for T = 4, and 
more data will be needed for SU(3) case. 



Figure 2. Error as a function of number of MC 
sweeps at T = 2 for U(l) = 1.2, SU(2) = 3.0, 
SU(3) = 6.25 and improved action for SU(3) 
= 3.3. 



3. Transport Coefficients 

We determine the parameters in the spectral 
function by making the fit of Matsubara Green's 
function. The fit of Gu(T) is shown in Fig. 3 for 
SU(2) at = 3.0, where the parameters and their 
errors are determined by the jackknife method 
using the least square package SALS. In this work 
we excludes the points where Error/ < G(T) > 
larger than unity from the fit range, and the data 
at T = is not included. 

We could not find the difference between the 
definitions of field strength tensor from the pla- 
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quette variables('diagonal' and 'perturbative' def- 
initions) in the present accuracy of data. Then at 
this (3 region, the higher order effects in ga 2 may 
be small for SU(2) gauge theory. 



Figure 3. Fit of Gn(T) by the parameters of 
spectral function for SU(2) at (3 = 3.0 



In the rcgin T = 3 — 4, Gu(T) becomes very 
flat, which is also observed for U(l) and SU(3). 
This behavior could not be found if we do the 
simulation in the smaller lattice of Nt — 4. 

The sheer viscosity rj and bulk viscosity £ are 
separated by taking the following combination, 

7] = a(A,7,m)i 2 , 
4 

£ = a(A,7,m)i2 - -a(A,7,m)n, 

where the subscript means that the transport co- 
efficients determined from G^. The sheer and 
bulk viscosities are shown in Fig. 4 for SU(2) case. 
The errors are estimated by jackknife method. 

The calculation for the SU(3) case is now un- 
der simulation at (3 = 3.3 and 3.2. The re- 
sult for the transport coefficients at f3 = 3.3 
are r\ x a 3 = 0.93 x 10~ 3 ± 0.107 x 10~ 2 , and 
£ x a 3 = 0.10 x 10~ 2 ± 0.13 x 10~ 2 . The results 
have still large errors and are very preliminary, 
but we think that within a few month the result 
will be improved and the data at (3 = 3.2 will also 
be presented. 



Figure 4. The sheer and bulk viscosity of SU(2) 
case as a function of (3 



We have found that, in the confined phase, it is 
very difficult to reduce the error of G^. We in- 
terpret it because the energy momentum tensor 
in the confined phase should be written by the 
hadron fields. In the deconfined phase the fluc- 
tuation for G^ v is still large, but by the use of 
Iwasaki's improved action, it is much suppressed 
and the transport coefficients for SU(3) gauge 
theory is calculated on 16 3 x 8 lattice. How- 
ever there are many thing to do before we obtain 
the result for the transport coefficients which are 
used for the quantitative phenomenology of quark 
gluon plasma. 
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Fig. 1: 

Error as a function od number of MC sweeps in the 
case of U(1) for confined and deconf iend phase. 
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Fig.2: 

Error as a function of number of MC sweeps at T=2 
for U(1) 3=1.2, SU(2) 3=3.0, SU(3) 3=6.25 and 
improved action for SU(3) 3=3.3. 
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Fig. 3: 

Fit of G1 1 (T) by the parameters of spectral 
function for SU(2) at 3=3.0. 
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Fig.4: 



The Sheer and Bulk Viscosity of SU(2) gauge theory 
as a f unct i on of 3 . 



